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ABSTRACT 



O ■ 

<n : 

5h ■ 

^ 1 We present the Spitzer/ Infrared Spectrograph spectrum of the main-sequence 

^ ■ star HD165014, which is a warm (> 200 K) debris disk candidate discovered by 

the AKARI All-Sky Survey. The star possesses extremely large excess emission 

PLh ■ at wavelengths longer than 5 /im. The detected flux densities at 10 and 20 /im 

are ~ 10 and ~ 30 times larger than the predicted photospheric emission, respec- 
tively. The excess emission is attributable to the presence of circumstellar warm 
dust. The dust temperature is estimated as 300-750 K, corresponding to the dis- 
tance of 0.7-4.4 AU from the central star. Significant fine-structured features are 
seen in the spectrum and the peak positions are in good agreement with those of 
crystalline enstatite. Features of crystalline forsterite are not significantly seen, 
i ■ HD 165014 is the first debris disk sample that has enstatite as a dominant form 

of crystalline silicate rather than forsterite. Possible formation of enstatite dust 
from differentiated parent bodies is suggested according to the solar system ana- 
log. The detection of an enstatite-rich debris disk in the current study suggests 
■ the presence of large bodies and a variety of silicate dust processing in warm 

q ■ debris disks. 

Subject headings: circumstellar matter — zodiacal dust — infrared: stars - 
stars: individual (HD 1650 14) 
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Introduction 



Debris disks were discovered in main- sequence stars by infrared excess over the pho- 
tospheric emission in observations by IRAS in the 1980s. Since debris disks are thought 
to be "extra-solar zodiacal light" formed via dust production from asteroids and/or comets 



(e.g. iBackman fc Parescd Il993t iLecavelier Pes Etangs et al.lll996l ). it is interesting to ex- 
amine mineralogical characteristics of debris dust and to explore connection between the 
debris dust and the dust in the solar system. Debris disks are also important as a probe of 
planetesimals, building blocks of planets, in extra-solar systems. 

Mid-infrared (MIR) spectroscopy is a strong tool to investigate th e properties of de- 
bris dust. Recent MIR observations w ith InfraRed Spectrograph (IRS; iHouck et al.l 120041 ) 
on board Spitzer ( IWerner et al.l 12004 ) have revealed the pres ence of abundant crystalline 



forsterite (Mg 2 Si04) and silica (SiC^) 
2009h . 



in several debris disks (e.g. lChen et al.ll2006l ; iLisse et al. 



In this Letter, we present an MIR low resolution spectrum of the main-sequence star 
HD165014 obtained with Spitzer/IRS. HP 16 5 14 is a debris disk c andidate with large MIR 
excess selected from a search for deb ris disks (IFuiiwara et al.ll2010l ) in the ^4/C4i?//Infrared 
Camera (IRC) All-Sky Survey data (jlshihara et al.ll2010l ). We report the detection of abun- 
dant crystalline enstatite (MgSi0 3 ) dust compared to forsterite toward the star and discuss 
the origin of the debris dust around HD165014. 



2. Observations and Data Analysis 
2.1. AKARI/IRC all-sky survey 



The S9W (9 u rn) and L18W (18 urn) images of HP 165 14 were taken with the IRC 

as pa rt of the All-Sky Survey 



( IQnaka et al.l 120071 ) on board AKARI ([Murakami et al. 



observations from 2006 May to 2007 August (jlshihara et al. 



2007 



2010|). The 5a sensitivity for 



a point source per scan is estimated to be 50 mJy in the S9W band and 90 mJy in the 
L18W band, and the typical absolute uncertainty in flux density is ~ 3% for the S9W band 
and ~ 4% for the L18W band at present. HD165014 was observed twice in the survey with 
the interval of 6 months and the fluxes at the two periods agree with each other within 
the uncertainty, indicating no significant variations in the flux. The position accuracy and 
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spatial resolution is estimated to be better than 2" and 10", respectively, securely concluding 
that the AKARI source is associated with HD165014. HD165014 was selected as a candidate 
of wa rm debris disk with la rge 18 /im excess by a search in the AKARI/IKC All-Sky Survey 
data (IFujiwara et al.ll2010l ). 



2.2. Spitzer/lUS Observations 

IRS observations of HD165014 were made on 2008 November 5 (AOR ID 26122752). All 
the four low-resolution modules, Short-Low 2 (SL2; 5.2-7.7 /mi) and Short-Low 1 (SL1; 7.4- 
14.5 /xm), Long-Low 2 (LL2; 14.0-21.3 /mi) and Long-Low 1 (LL1; 19.5-38.0 /mi), were used to 
obtain a full 5—35 /im low-resolution (A/AA ~ 100) spectrum. We use the pipeline-processed 
(S18.1) Basic Calibrated Data products from the Spitzer Science Center and analyzed the 
sky-subtracted extracted one-dimensional spectral data "bksub.tbl" since the target is a point 
source in a relatively empty field. The wavelength calibration is as good as 0.1 /im along 
the dispersal direction. For the LL1 spectra, we use the data only for A < 35 /xm because 
the noise becomes large at A > 35 /im, whose spectral range is not crucial for the present 
analysis. The absolute flux accuracy in the SL and LL spectra of the pipeline-processed 
products is better than 10%. 



2.3. Subaru/COMICS Observations 



HD165014 was observed with the COoled Mid-Infrared Camera and Spectrometer (COMICS; 
Kataza et al.ll2000l ) mounted on the 8 m Subaru Telescope on 2008 July 16 and 17. Imaging 
observations in the 8.8 /im, 11.7 /im, and 18.8 /im bands were carried out. The secondary 
mirror chopping metho d was used for the background cancellation. We used standard stars 
(7 Aql and e Sco) from I Cohen et al.l (119991 ) as a flux calibrator and the reference point-spread 
functions were derived from the observations. We observed the standard stars before or after 
the observations of HD165014 in the same manner as HD165014. For the data reduction, 



we used our own reduction tools and IRAF (ITodyl 119931 ) . The standard chop pair subtrac- 
tion and the shift- and- ad d method in units of 0.1 pixel were employed. We applied airmass 
correction using ATRAN (lLordlll992l ). The derived flux densities are summarized in Tabled) 
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2.4. Photospheric Fitting 



The photospheric flux density of HD 1650 14 is estimated f rom the Kurucz model (IKurucz 



19921) fitted to the BVI- and Two Micron All Sky Survey flSkrutskie et all 120061 ) JHK S - 



band photometry of the star taking account of the extinction. The only a vailable spectral 



class ification of HD165014 (F2V) is provided by the Michigan Catalog (IHouk & Cowley 



19751 ). We reexamine the spectral type of the star by the photospheric fitting. We use 
the Kurucz models of various spectral-type dwarfs for the stellar photospheric template of 
HD165014 and the dis tance d and the extinction Ay are set free. We adopt an extinction 



curve A\/A v given by iFitzpatrick fc Massal (120091 ). 



A x _ / 0.349 + 2.087R V \ 1 

A v Vi + WO-SO?) 2 - 05 V )Ry U 

where A\ is the extinction toward the star at a wavelength A (in /im) and the ratio of the 
total to the selective extinction Ry is assumed as 3.11. 

The photospheric fitting shows that the B VI JH i^g-band photometry of HD165014 
is well accounted for with almost the same significance by spectral types between B8V 
(T eff = 12000 K) with Ay = 3.0 at d = 190 pc and F2V (T cS = 7000 K) with Ay = 1.9 
at d = 70 pc. It is difficult to determine the spectral type of the star more accurately only 
from photometric data due to its large extinction and optical high-resolution spectroscopy 
is required for further examination. An earlier spectral type at a larger distance seems more 
reasonable for HD165014, taking account of the large amount of extinction. In the following 
we assume A0V (with Ay = 2.7 and d = 140 pc) as the spectral type of HD165014 rather 
than F2V given by the Michigan Catalog. The estimated extinction Ay = 2.7 is still large 
for a star at 140 pc, suggesting that part of extinction may be of circumstellar origin. The 
derived photospheric flux density from the fitting in each MIR band, which does not depend 
on the spectral type sensitively (only a few percent of variation from B8V to F2V), is listed 
in Table [TJ 



3. Results 

3.1. Spectral Energy Distribution 

The obtained Spitzer/IKS spectrum of the star is shown in Figure [1] together with the 
AKARI and Subaru/COMICS photometry. The measurements of the MSX counterpart 
(MSX6C G009.0807+00.3009) and the Spitzer/IRAC GLIMPSE II counterpart (SSTGLMC 
G009. 0801+00. 3011) are also shown in Figured! Most of the MIR flux densities are in good 
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agreement with each other, while those of AKARI 18 /im band and Spitzer/IKAC 8 /im-band 
are slightly higher than the others. The expected photospheric emission is also plotted in 
Figure [TJ Significant excess emission at wavelengths longer than 5 /im is clearly seen. The 
detected flux densities at 9 and 18 /im are 10 and 30 times larger than the photosphere, 
respectively, unambiguously indicating the presence of a warm and bright debris disk. The 
slope of the observed spectrum at A > 20 /im is consistent with that of the Rayleigh- Jeans 
side of a blackbody, suggesting the dust temperature of > 300 K and the truncation of the 
disk at the radius where the dust temperature is ~ 300 K. 

One of the observable indicators of dust abundance is its fractional luminosity, a ratio of 
the infrared luminosity from the disk to the stellar luminosity. The fractional luminosity of 
HD1 65014 is estimated as ~ 5 x 10~ 3 . It is comparable to (5 Pic fearrado y Navascues et ah 



19991 ). suggesting that HD165014 is one of the brightest debris disks discovered so far. The 
location of the star, deep in the galactic plane, makes it difficult to find until the AKARI 
survey. 



3.2. Features in the Spectrum of the Excess 



To examine the excess emission, we subtract the expected photospheric emission from 
the observed IRS spectrum. The photosphere-subtracted spectrum shows significant emission 
features centered at around 10 and 20 /im. A ubiquitous dust species family, silicate, which 
has features around ~ 10 /im (due to Si-0 stretching modes) and ~ 20 /im (due to O-Si- 
O bending modes), seems to be a main carrier of the observed features. Particularly, the 
spectrum shows several fine-structured features attributable to crystalline silicates. 

To investigate the details of the features, Figure |2] plots the normalized emissivity of 
the excess emission of HD165014 by dividing the photosphere-subtracted spectrum by a 
blackbody (B V (T)) of a temperature T, where T is chosen so that the spectrum level of the 
feature-free region (8.0 and 13.2 /im for the top panel, and 13.5 and 33 /an for the bottom 
panel of Figure[2]) becomes almost flat. The chosen T is 560 and 330 K for the top and bottom 
panels of Figure El re spectively. For co mparison, spectra of crystalline enstatite, forsterite 
(ITamanai et al.l 120061 ) , and fused silica (IKoike et al.lll9891 ) are also plotted in Figure [2j 



A significant broad feature is seen at 9-12 /im in the spectrum. Since the width of this 
feature is large (FWHM ~ 2.8 /x m), it seems to originate from /zm-sized amorphous silicate 
rather than sub-/im-sized silicate (Ivan Boekel et al.l 120031 ) . In addition, some peaks are also 
seen on the broad feature. The most significant peak is located at 9.3 /im, and the second at 
10.5 /im. Weaker narrow features are also seen around 11.1-11.5 /im. The positions of all the 
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narrow features are consistent with the spectrum of crystalline enstatite. The 9.8 /mi feature 
of crystalline forsterite, which is the most common form of crystalline silicate, is not seen in 
the excess emission of HD165014. Only weak narrow features are seen around 11.1-11.5 /mi, 
where a strong crystalline forsterite feature is expected, in the spectrum of HD165014. 

To analyze the iV-band spectrum quantitatively and confirm the identification of the 
features described above, we make a simple spectral model fitting of the deriv ed emissivity at 



8-13 /mi. We consider four dust species, cry stalline forsterite and e nstatite (ITamanai et al. 



20061 ) and amorphous olivine and pyroxene (IDorschner et al.lll995l ). in our model fit. The 



model spectrum is given by 

Emissivity = a cont + a iQ i + a px Q px + a io Q io + a cn Q en , (2) 

where Q \, Q px , Qf a , and Q cn are the absorption efficiencies of amorphous olivine, amorphous 
pyroxene, crystalline forsterite, and crystalline enstatite, respectively. The scaling factors, 
as, are free parameters (a con t is the contribution of continuum). For amorphous olivine and 
pyroxene, we choose one dust size from 0.1, 1. 0, 1.5, or 2.0 /tm , whose absorption coefficients 



are computed from the op tical constants of IDorschner et al.l (119951 ) based the Mie theory 



(IBohren fc Huffman! 1 19831 ). We employ a least-squares minimization method to determine 
the most likely parameters. A model with 2.0 /tm sized amorphous silicates is found to 
provide the best fit among all the considered dust sizes of amorphous silicates. The derived 
best-fit parameters are a cont = 0.44, a i = 0.00, a px = 0.24, af = 0.03, and a en = 0.32. The 
best-fit model spectrum is shown in the top panel of Figure [31 indicating that the iV-band 
spectrum of HD 1650 14 is well reproduced by the model spectrum. The broad and narrow 
features seen in the observed spectrum are indeed attributed to 2 /tm sized amorphous 
pyroxene and crystalline enstatite, respectively. The contribution of forsterite dust is small 
and the mass ratio of forsterite to enstatite is estimated as ~ 1/20. 

In the longer wavelength region of the spectrum, a significant trapezoidal feature is seen 
at 17-20 /mi, which is very similar to the enstatite spectrum. A broad feature seen around 
28 fim also seems to originate from enstatite. A weak feature is seen around 23 /tm, which 
may be accounted for by both forsterite and enstatite. The feature around 23 /tm might be 
attributable to forsterite; forsterite, if any, is less abundant than enstatite since the feature 
is very weak. Silica features at 9.0 and 20-21 /tm are not seen in the excess emission of 
HD165014. 

The observed features and their identifications are summarized in Table [2J From the 
overall shape of the MIR excess spectrum of HD165014, we conclude that most of the fine 
features are attributable to crystalline enstatite and that forsterite is much less abundant 
than enstatite around the star. The feature-to-continuum ratios in the 10 and 20 /mi regions 
are ~ 1 and ~ 0.25, respectively. The contribution from the continuum component is large 
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in the excess emission of HD165014, suggesting the presence of abundant featureless large 
dust grains (> 10 /im in size) around the star. 



3.3. Radial Distribution of Dust 

In the top panel of Figure [21 a rising continuum toward shorter wavelengths A < 7 /im 
is seen, suggesting the presence of hotter (T > 560 K) materials around the star. Although 
the rising continuum might also be attributable to amorphous carbo n, whose emissivity rises 



toward shorter wavelengths in the NIR and MIR flJager et al.lll998[ ). no more discussion can 
be given with the present data. The temperature of the hottest component is estimated 
as T = 750 K from a fitting of the 5 — 7 /im excess spectrum with a single-temperature 
blackbody assuming that the 5 — 7 /im excess is coming from the innermost hottest region. 
Therefore dust grains are distributed in the region with temperatures of T = 750 — 300 K 
around HD165014. The inner and outer radii of the debris disk are estimated as 0.70 AU 
and 4.4 AU, respectively, by assuming that the debris dust is a blackbody particle. Possible 
large circumstellar extinction suggests that the disk is seen nearly edge-on. 



4. Discussion 



As mentioned above, it is evident that crystalline enstatite is abundant in the d ebris disk 
around HD 1650 14 and forsterite features are not clearly seen. IChen et al.l ( 120061 ) conduct 
a comprehensive MIR spectroscopic survey of debris disks. Among their 59 debris disk 
samples, five objects (HR3927, rj Crv, HD1 13766, HR7012, and rj Tel) are found to possess 
spectral features that are well-modeled by /zm-sized amorphous and crystalline silicates. 
Forsterite features are dominant compared to enstatite toward four samples (HR3927, rj Crv, 
HD113766, and rj Tel). HR7012, which shows a very stro ng peak at 9.1-9. 2 /im in its 
spectrum, was in itially modeled by enstatite- rich silicate by IChen et al.l ( 120061 ). However, 
Lisse et al.l (120091 ) conclude that the observed 9.1-9.2 /im feature originates from abundant 
silica grains and that the feature strengths of enstatite and forsterite are almost comparable. 
Although there are a few additiona l debris disk stars with significant dust features in the N - 
band (p Pic, HIP8920, HD145263; iKnacke et al.lll993[ ISong et~al]l2005[ iHonda et alJliooi ). 
no sample with high abundance of enstatite is reported to date. HD165014 is the only known 
debris disk that shows strong enstatite features compared to forsterite. 



In the solar system, a number of achondrite meteorites, which are primarily composed 
of enstatite, named Aubrites, have been discovered. E-type asteroids are thought to have en- 
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statite achondrite surfaces and to be parent bodies of aubrites (e.g. IZellner et al.lll977l ) based 
on their reflection spectra. E-type asteroids form a large proportion of asteroids inward of 
the main belt known as Hungaria asteroids and are bel ieved to be frag ments of differenti- 
ated larger bodies that were heated at least to 1700 K ( IKeil et al.lll989[ ). Considering that 
enstatite is dominant around HD 1650 14 rather than forsterite, the observed debris dust may 
originate from an analog of E-type asteroids in the solar system, being h armonic with the idea 
that the debris disk is formed by means of collisions between aster oids (IBackman fc Paresce 
1993 ). It is also known that Mercury's surface is rich in enstatite ( Sprague fc RoushT 19981 ). 
If a Mercury-like planet exists around HD 1650 14 and a mechanism to scatter the surface 
material of the planet, for example infall of small bodies, occurs, an enstatite- rich debris 
disk might be formed. It should also be noted that Mercury's high de nsity is interpreted 
as a result of stripping of the surface crustal material (IBenz et al.l ll988). The possible link 
between the debris material around HD165014 and Mercury needs to be further explored. 

Although a number of MIR spectra of protoplanetary disks associated with younger 
stars, Herbig Ae/Be and T Tau stars, have been collected so far, objects that show enstatite 
features are rare. HD179 218 is a Herbig Ae/Be star, which has the most enstatite-rich 



spectrum known to date (jBouwman et al.l 12001 



Schiitz et al 



2005J)- W e show the iV-band 



spectrum of HD179218 obtained with TIMMI2 (Ivan Boekel et al.ll2005f ) in the bottom panel 
of Figure [3j While crystalline enstatite features are clearly seen in the spectrum, features 
at 9.85 and 11.20 /zm attributable to forsterite are also seen with comparable strengths 
as enstatite, suggesting the presence of abundant crystalline forsterite as well as ens tatite. 
Annealing experiments of a magnesium silicate smoke made by iRietmeijer et al.l (119861 ) show 
that the initially formed forsterite and silica react with each other and form enstatite by the 
following reaction 



Mg 2 Si0 4 + Si0 2 — ► 2MgSiO 



3- 



(3) 



Bouwman et al.l (120011 ) suggest that the presence of enstatite around HD179218 might be 
due to the high luminosity (3OQL ) , which gives rise to rap id dispersal of the gas, resulting 
in a high efficiency of the reaction. Ivan Boekel et al.l ( 20051 ) suggest that enstatite might be 
produced by means of chemical equilibri um processes in h igh temperature environments, i.e., 
inner regions of the protoplanetary disk. ISato et al.l (120061 ) show that crystalline enstatite can 
be produced by simultaneous evaporation of SiO grains and Mg vapor in a plasma in labora- 
tory experiments. It is now widely a ccepted that the turbul ence in disks is attributed to the 
magnetorotational instability fMR I; iBalbus fc Hawlev!ll9911) . suggesting that at least part of 
the disk is sufficiently ionized (e.g. Inutsuka fc Sand 20051 ). Therefore the formation process 
of crystalline enstatite suggested by lSato et al. ( 2006 ) may work efficiently in protoplanetary 
disks. It is possible that crystalline enstatite grains produced in the protoplanetary disk were 
once stored in small bodies such as comets and have been released recently in the HD 165014 
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system. Indeed crystalline enstatite is detected towar d some comets in the solar system by 



astronomical MIR o bservations (e.g. iLisse et al.ll2006l ) as well as by the STARDUST sample 



return mission (e.g. IZolensky et al.l 120061 ). However, it remains an open question why crys- 
talline forsterite is depleted around HD165014. Further experimental and theoretical studies 
are required to discuss the origin of abundant crystalline enstatite in the debris disk around 
HD165014 in detail. 
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Table 1: Infrared Photometry of HD165014. 



A 




Instrument 


Photosphere 11 


Significance, x b 


w 


(Jy) 




(Jy) 




5.8 


0.90 ±0.02 


Spitzer/IRAC C 


0.26 


27 


8.0 


1.06 ±0.02 


Spitzer/IRAC C 


0.14 


43 


8.28 


1.08 ±0.05 


MSX d 


0.14 


22 


8.8 


0.93 ±0.09 


COMICS 


0.12 


10 


9 


1.36 ±0.10 


AKARI/IRC 


0.14 


13 


11.7 


0.89 ±0.09 


COMICS 


0.069 


10 


12.13 


0.78 ±0.07 


MSX d 


0.064 


11 


14.65 


0.76 ±0.06 


MSX d 


0.044 


13 


18 


0.96 ±0.14 


AKARI/IRC 


0.030 


7 


18.8 


0.84 ±0.01 


COMICS 


0.027 


220 



"From the Kurucz model to fitted to BVIJHK s -b&nds data. 
b X = (Observed — Kurucz)/noise. 

c Thc Spitzer GLIMPSE II counterpart (SSTGLMC G009.0801+00.3011), whose positional offset from 
HD165014 is O'.'IO. See http://www.astro.wisc.edu/glimpse/glimpsedata.html, 
d The MSX counterpart (MSX6C G009.0807+00.3009), whose positional offset from HD165014 is 2'/ 58. 
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Table 2: Features in Excess Emission of HD165014. 



A 


Comment 


Identification 


(/im) 






10 (7.3-12.4) 


Very strong broad feature 


Silicate 


9.3 


Strong peak within 10 /im feature 


Crystalline En 


10.5 


Strong peak within 10 /im feature 


Crystalline En 


11.1 


Peak within 10 /im feature 


Crystalline En/Fo 


11.5 


Shoulder within 10 /im feature 


Crystalline En 


14.0-15.0 


Weak feature 


Unidentified 


18 (17.8-19.7) 


Strong trapezoidal emission 


Crystalline En 


22.8-24.3 


Weak feature 


Crystalline En/Fo 


27.3-30.3 


Weak feature 


Crystalline En 


33.8-34.5 


Possible weak feature, low signal-to-noise ratio 


Unidentified 
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Fig. 1.— Spectral energy distribution of HD165014. The open squares, triangles, cir- 
cles, crosses, and pluses indicate the photometric data obtained with AKARI/IRC, Sub- 
aru/COMICS, MSX, Spitzer/IRAC, and BVIJHK s -ba,nd photometry taken from Vizier 
database, respectively The solid and dotted lines indicate the Spitzer/IRS spectrum and 
the photospheric contribution of an A0V star (jKuruczlll992l ). respectively. 
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Fig. 2. — Top: emissivity of the excess emission of HD165014 at 5-15 /xm derived by division 
of the observed spectrum by a blackbody of T = 560 K (thick solid line). For comparison, 
spectra of crystalline enstatite, forsterite, and fused silica from laboratory measurements are 
also plotted as thin solid, dashed, and dotted lines, respectively. Bottom: emissivity of the 
excess emission of HD165014 at 12-35 /im derived by division of the observed spectrum by 
a blackbody of T = 330 K. The line styles are the same as the top panel. 
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Fig. 3. — Top: fit result of the iV-band emissivity of HD165014. The thick solid line is the 
best-fit model spectrum, which is the sum of constant continuum (thin solid line), 2.0 /jm- 
sized amorphous pyroxene (dashed line), crystalline forsterite (dotted line), and enstatite 
(dot-dashed line). Amorphous olivine component is not shown since no contribution of it 
is required in the best-fit model. The residual spectrum subtracted by the best-fit model is 
shown below the fit result. Bottom: comparison of the iV-band spectra of HP 1650 14 (thick 



solid line) and enstatite-rich Herbig Ae/Be star HD179218 (thin solid line; Ivan Boekel et al. 



20051 ). The vertical dashed and dotted lines indicate peak positions of crystalline enstatite 
(En) and forsterite (Fo) detected toward HD179218, respectively. 



